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Theoretical predictions for Bhabha scattering at the two-loop level require the inclusion of 

f^ ■ hadronic vacuum polarization in the photon propagator. We present predictions for the con- 

00 ! tributions from reducible amplitudes which arc proportional to the vacuum polarization 7r(g^) 

^D ' and from irreducible ones where the vacuum polarization appears in a loop representing vertex 

or box diagrams. The second case can be treated by using dispersion relations with a weight 

function proportional to the i?-ratio as measured in electron-positron annihilation into hadrons 

. , and kernels that can be calculated perturbatively. We present simple analytical forms for the 

j^ ' kernels and, using two convenient parametrizations for the function R{s), numerical results for 

the quantities of interest. As a cross check we evaluate the corresponding corrections resulting 

from light and heavy lepton loops and we find perfect agreement with previous calculations. For 

the hadronic correction our result are in good agreement with a previous evaluation. 
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1 Introduction 

Electron-positron colliders, with their potential for precise and specific measurements of cross sections, 
have been and are being operated from the very low energy region around the pion threshold up to more 
than 200 GeV and may, in the future, reach up to energies of one or perhaps even several TeV. To determine 
the luminosity, one necessarily uses a reaction, whose cross section can be well measured and furthermore 
calculated with sufficient precision. Reactions which involve only leptons in the final state, like electron- 
positron annihilation into muon pairs or elastic electron-positron (Bhabha-)scattering are ideally suited for 
this purpose. In particular Bhabha scattering, with its relatively large cross section, has always been the 
standard luminosity monitor reaction. Precise theory predictions are, therefore, mandatory and, in view of 
recent interest in precise measurements with high counting rates, must be pushed to two-loop order. 

Various ingredients are necessary in this connection. A major step has been made in [1] where the 
photonic two-loop virtual corrections plus the corresponding soft real radiation has been evaluated for the 
case of interest to^/s ^ 1, using earlier results for the completely massless case [2] and exploiting the 
relation between the soft and coUinear singularities for these two limiting cases. These results were confirmed 
in [3], where in addition also the contributions from muon loops (again in the limit m? / s <C 1 ) were 
calculated. These muon-loop contributions, in the same high-energy limit were also evaluated in [4j, the 
corresponding results for arbitrary mass of the internal lepton were presented in [5]. The electron loop 
corrections involving the exact m^/s dependence were computed in [6], while further efforts towards the full 
electron mass depedence at two-loop level can be found in 7 . 

All these contributions can be calculated strictly within Quantum Electrodynamics (we do not consider 
electroweak corrections that are relevant for high energies). However, in two loop approximation contribu- 
tions from virtual hadrons come into play. From general considerations it is obvious that, generally speaking, 
their magnitude is comparable to or larger than those from virtual muons. It is well known jSj that these 
virtual hadronic contributions can be evaluated through dispersion relations, folding the absorptive part 
of the hadronic vacuum polarization, the i?-ratio measured in electron-positron annihilation, with a kernel 
that can be calculated perturbatively, in the present case in a one-loop calculation. This approach has been 
adopted in [9], where the kernel for the box diagram has been calculated for non- vanishing electron mass 
and the limit jrie -^ has been considered only subsequently. A more compact form for this kernel can be 
obtained by using directly the well known results for the (direct plus crossed) box with one photon and one 
massive vector boson, the Z-7 box, contributing to Bhabha scattering at high energies. Since it is well known 
that in this case the electron mass can be safely set to zero from the beginning , the calculation becomes 
significantly simpler. Vertex corrections with a hadronic insertion have been evaluated with this technique 
long time ago and it is only this box contribution, that was not yet available since long. 

As stated above, results for the hadronic contributions have been presented in [9l. In view of the fact, 
that we are using a somewhat different approach and furthermore, to provide an independent cross check 
the results of our calculation will be presented in some detail. To allow for an easier comparison with [^, 
we shall use the same parametrization [lOj for the i?-ratio. In addition we shall compare the results to the 
ones derived from a second parametrization [11| that includes more recent data and these latter ones should 
be considered as our definite predictions. 

The paper will be organized as follows: In chapter 2 we present the general analysis and classify the various 
reducible and irreducible contributions. In chapter 3 we give the details of the calculation, the explicit forms 
of the kernels and identify the contributions from real radiation needed to render the results infrared finite. 
In this connection it is convenient to split the virtual (plus soft real) corrections into different building blocks 
that will be described in more detail in this chapter. The handling of the dispersion integrals with their 
poles is described in chapter 4. Of specific interest is the high energy limit, with s, t and u in an region 
where the R- ratio has approached an approximately constant value. In this limit and in analogy with the 
treatment on the form factor in 12J a particularly simple form can be derived where the information about 
hadron physics can be encoded into three "moments" of the R ratio. Using this method and evaluating the 
moments for a lepton, e.g. the muon or the tau-lepton, the results from ^3j are easily recovered. In chapter 
5 we present the numerical results for the two parametrizations. We give the results for the building blocks, 
vacuum polarization, vertex and boxes, and the complete corrections, split into the various contributions 



and discuss their physical relevance. The final Chapter 6 contains a brief summary and our conclusions. 

2 General analysis 

It is well known that contribution from the hadronic vacuum polarization can be directly evaluated by 
convoluting an appropriately chosen kernel with the familiar _R-ratio (R = (Thad/cpt) measured in electron- 
positron annihilation [8j. An arbitrary amplitude involving the hadronic vacuum polarization is obtained, 
by definition, from the original one by replacing the photon propagator as follows: 

2 , ^ -T—-nq9 - g g ) n(g ) ': . (i) 

The renormalized vacuum polarization function 11 (g^) is obtained from its absorptive part (essentially the 
i?-ratio) by the subtracted dispersion relation: 

n?" =-- / 2 T^' ImUiz) = --R{z), 2 

"■ J4m2 z q-' - z + le 3 

and has a cut for g^ > 4m^, with the threshold for hadron production at 4to^. The q^q'^ term in Eq.(IT]) does 
not contribute and the photon propagator is effectively replaced as follows: 

-igafj -igaf}^, 2n ■ a f°° dz R(z) 

q'^ + le g^ + le 6tt J^rn^ z q^ ~ z + le 

If q^ is fixed by the external kinematics, applying the correction is equivalent to multiplication of the previous 
amplitude by Il{q^). In higher orders, summing the one-particle reducible terms only, this corresponds to 
the replacement of the photon propagator by the dressed one: 



Wa/3 - igap 1 
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(4) 



However, if q stands for a loop momentum, it is convenient to exchange the order of integration and evaluate in 
a first step the loop integral with a ficticious massive vector boson of mass y^, and to convolute subsequently 
this amplitude with the i?-ratio, i.e. with -§;;!^2^ R{z) ■ In [12 this has been done for the Dirac form factor, 
assuming massless external fermions, and special emphasis has been put on the investigation of the limit, 
where the momentum transfer is far larger than A:m? . A similar approach will be useful for the present case. 
We will, in a first step, investigate the generic case valid for arbitrary s, |t|, |u| ^ Wg. The weight functions 
needed to obtain the three building blocks, namely the vacuum polarization function n(g^) (Fig. [TJi), the 
correction to the vertex function V{q^) {Fig-Wp) and the amplitude arising from the box diagrams (Fig. [TJ;), 
can be taken from the literature. 
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Figure 1: The building blocks for the QED hadronic corrections 

The complete hadron induced corrections are conveniently split into thre classes: 

1) Tree level diagrams, with two vacuum polarization insertions proportional to n(s) or n(i) where 11 
originates from virtual hadrons, muons or electrons (Fig. [2]). 

These corrections to the amplitudes are proportional to n(s)^, n(t)^ or n(s)n(t) (reducible quadratic 
Il-terms) and are directly obtained from the Born amplitude. 




2 Re 




>K"X"^^^ 




Figure 2: Tree level diagrams with vacuum polarization insertion 
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Figure 3: One- loop photonic diagrams in combination with a dressed photon propagator 

2) Corrections which involve one-loop purely photonic corrections in combination with a dressed photon 
propagator in the s or t channel (Fig.[21). 

These are proportional to n(s) or n(t) (reducible linear n-terms) and are directly obtained from the 
one-loop corrections to the Bhabha scattering. They can again be separated into amplitudes resulting 
from one- loop vertex and box corrections, respectively. Both are rendered infr'ared finite by adding soft 
real photon emission with E-.^ < lu << ^/s. In both cases this leads to a logarithmic a;-dependence. 

The photon vertex correction involves a coUinear electron-mass singularity, which leads to the only rrie 
dependence relevant for our investigation. 

Also the 77-box amplitude, after interference with s- or i-channel dressed photon exchange, leads to a 
logarfthmic dependence on oj. The electron mass, however, may safely be set to zero. 

3) As a third class we have to consider the irreducible two-loop contributions, i.e. amplitudes with dressed 
photon propagator in a loop (Fig.|4]). 
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Figure 4: Irreducible two-loop diagrams 



The dressed propagator may be located either in a vertex or in a box, which wiU interfere with the 
two Born amphtudes from s- and i-channel exchange. Irreducible vertex corrections are infrared finite 
and me may be safely set to zero. They are easily obtained from the Born amplitude by replacing the 
appropriate vertex by V{s) or V{t) as defined below. The box amplitudes are again infrared divergent 
and must be made finite by combining with soft real radiation. 

Let us emphasize that contributions from lepton loops follow the same classificatiorlj. The simple form of 
the result in the high energy limit, s, |i|, |u| ::^ 4m J, will allow for a convenient cross check of our calculation. 
In addition to these corrections with virtual photons, one has to compute the corresponding emission of 
real photons (Fig. [5|) to compensate the infrared divergencies. 
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Figure 5: Diagrams with real photon emission. The four- momentum of the real soft photon is fc = (fco,fc) and 
the integration is performed for \k\ < to. Here and below infrared regularization through a small photon mass A 
(A^ = ko — k^) is implicitly understood. 



3 Details of the computation 

In this section we discuss these contributions in details. In all equations containing products of Feynman 
diagrams, the sum over the spins of the outgoing particles and the average over the spins of the incoming 
particles is implicit, as well as the conservation of the external momenta. In these formulas, the coefhcient 
Cps = (64s7r^)~^ comes from the integration of the phase-space of the outgoing electron and positron. 

3.1 Vacuum polarization insertion 

The Born cross section is obtained from the combination of s- and i-channel exchange: 
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Replacing the photon propagator in the s- and i-channel by the dressed one, one obtains: 
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(6) 



with X — —t/s = (1 — cos6')/2 and 



n(g2) = ne(g2) + n^(g2) + n,((72) + nhad(g'). 



(7) 



^With the exception of the two-loop vacuum polarization (Fig.[6ll whose absorptive part is in the hadronic case, by definition, 
part of the H-ratio. 



Expanding up to order o? ^ one easily recovers the Born contribution and the reducible one- and two-loop 
corrections: 
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Using for light leptons (e and ^) 
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the well known electron/muon induced one- and two-loop reducible contributions are easily recovered. In the 
present context the two-loop terms involving hadrons arise from terms proportional to H^^^^ and nhadlliopt, 
with different combinations of real and imaginary parts. 

3.2 Reducible diagrams 

Contributions from one-loop photonic amplitudes, interfering with amplitudes with the dressed photon 
propagator in the s- or t-channel are infrared divergent and must be combined with real radiation. For the 
amplitudes involving vertex corrections with have: 
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The sum of these three contributions from s-channel, t-channel and their interference gives the differential 
cross section: 
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where A is the photon mass used as IR regulator. We have also introduced 
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Just like for the one-loop corrections, a logarithmic dependence on me from coUinear singularities remains. 
It is interesting to notice the presence of an infrared divergent term proportional to Imn(s) surviving after 
the inclusion of the real soft photon emission. Remembering that these contributions can be easily obtained 
from the 0{a) calculation it is clear that the 0{a) result is easily recovered by the substitution of n(s) and 
n(i) by 1. 

A similar discussion applies to the photonic box diagrams, interfering with amplitudes with the dressed 
photon propagator in the s- or i-channel: 
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The differential cross section is then given by: 
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In this case the electron mass can be safely set to zero. Also in this case a logarithm of A^ survives and cancels 
exactly the one generated by the vertex reducible corrections, rendering the total cross section infrared finite. 



3.3 Irreducible diagrams 

Let us move to the third group consisting of the two irreducible contributions. The vertex correction has 
been discussed in detail in 12J and can be cast into the following form: 
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The contribution to the cross section can be cast into a form closely related to the Born cross section. 
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As stated above, m^ has been set to zero and the result is obviously infrared finite. 

Finally for the irreducible two-loop box contributions the kernels can again be directly taken from the 
literature '13] . The part of the kernel, which corresponds to the infrared divergent piece will be canceled by 
the proper combination of real soft radiation amplitudes which are also proportional to n(q^), specifically: 
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In total we find: 
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where 
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The part proportional to ln(2aj/-\/s) has been displayed separately and, as stated before, is proportional to 
n(s) or n(i). As discussed in the introduction, the functions S^a and ^B^ corresponding to the box diagram 
with a massive and a massless vector boson, can be directly read off from the literature |13j . 



For the two-loop hadronic contributions to cross section we find (without the trivial vacuum polarization) 
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4 Evaluation of the dispersion integrals 

In Eq. (|5D|) . the total cross section is written in terms of the building blocks Il{q^), ¥{0^) and B{a,b,c) 
defined in Eq.®, Eg. ipD]) and Eq. pS)) respectively. In the hadronic case, given a suitable parametrization 
of R{s), these dispersion integrals have to be integrated numerically. Therefore, before attempting the 
evaluation of the cross section, all sources of numerical instability must be cured. 

The expression for n(g^) in Eq.([2]) is very simple, but reveals the presence of a pole of the integrand 
at z — q^ + ie. The simplest way to get rid of it is to add and subtract R{q^) in the integrand for g^ > 
(s-channel). After the useful change of variable z = Am^/y, we get: 
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The integral for V{q'^) given in Eq. ipO)) does not show any pole in the integration domain and is directly 
accessible to a numerical evaluation. However, its convergence in the high energy integration region can be 
improved introducing the asymptotic, approximately constant value R{oo) of the i?-ratio. To this purpose, 
let us recall the results from '12] for the vertex V, which can be rewritten in the form: 
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The first one of these integrals can be solved exactly: 
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where r = q^ /(Am} — ie), and the second integral converges well in the large momentum region. 

A similar approach can be adopted for integrating the kernel from the box diagram B{a, b, c) defined in 
Eg. ([29| . The function £^a has a good high energy behaviour, but has a pole (for a = s > 0) in z — a + ie and 



can be treated in the same way as Il{q^): 
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The first integral in the expression for Ba{s, b, c) is then given by: 
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In the last expression we have introduced x = 1 — x, r = s/{Am? — ie), r = 1 — r. On the other hand, Bi 
can be computed following the same procedure used for V(q^): 
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where for the first integral we have: 

lo y y 



lB{a,b,c)= / — S,Bia,b,c, ) 

Jo y y 



Li2 ( ^ I — In T- 

Amr I Amr 



^ m--^ 



Am? 



\Am? 



yAm?) 



\Arr?J ' \An? 



Am? — a 



—b c\ An? — a t ■ ( '^ \ 
An? a J Am"^ \Am?J 



Am^ + b 



-b An? + b . ( -b\ 
An? Am"^ \An?J 



(38) 



4.1 High energy limit 

In the high energy limit, i.e. for ^/s and \/— t far larger than the energy above which R{s) approaches 
(sufficiently rapidly) i?(oo), the building blocks Hiq^), Viq^) and B(a,b,c) can be expressed in terms of the 
moments i?„ defined through: 
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The large q^ behaviour of the vacuum polarization is then given by 
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and V takes the following form: 
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-2^^ 



(41) 
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Similarly, the building block B{a, b, c), in the high energy limit s, |t|, |u| 3> '^m? is given by: 



B{a,h,c) = — { i?(oo) I^ {a,b,c) + Rq 



c? '"'[ Am? 2 



c-b 
"2 



-(4+6C(2))+-L6 
a \ /a 



Ri—Lbc (■ , 



(42) 



where 



I^[a,b,c) 



L, 



be 



C-b r 



2a 



1 2 ^"^ 1 ^'^ 

2 4to2 2 iw? 

a 



Li + 6C(2) 



In 



Arri? 



1 



(Lfc + Q + C(2) 



.7(5, c) + J(c, 6) 



4r?T? 



r2 



5C(2) 



+ 0(m2). (43) 



In the last equations we have introduced: 



Lh = In 



b+ie 



In 



c+ie 
a 



In- 



b+ie 



J{x,y) = 5"]; 



y-ie 



i nhi 



y-ie 



(44) 



4.2 The leptonic contribution 

With these ingredients the two-loop result induced by massive and light leptons is easily recovered. 

doj_^ d'^lii daf^d,i daj dcrf 

dfi dn dn dfi dn' 

All ingredients are obtained from the corresponding expressions for the hadronic case using the i?-ratio 



(45) 



Riiz) 



1 



Ami 

~2^ 



Amf 
1 ^. 



(46) 



Numerical evaluations for the muon and r-lepton will be presented below. In the high energy limit we will 
use the moments 1121: 



S 1 ^ 2K 

i?,(«^) = l, i?to-ln4--, i?z,i = -ln24--ln4+ — -C(2), 



Ri.2 = lln'A-lln\ + 2 
6 6 



|-C(2)]ln4 + 2C(3) + ^C(2)~^. 



(47) 



The integral over z can then be analytically evaluated, giving for the building blocks: 



nK9')--^( inT^-oK 



a ( 1 .-q^ 19, 2-9' , 1 



yii'i') = -iz i-j^'z^-^^^'z^ + -. 



B,(a,b,c) 



Svr I 12 
a 
3^ 



24 mf 



./„^ 265' 



, -r ..X 19 , , 3355 

1-7;^ +«^)- 12 ^(2) -13^ 



a 

c-b 
2a 



^bc 



In' 



In 



rn? 



Lh + Lr 



10 

y 



(ife + ic) 



56 



J(fe,c) + J{c,b) 



i6+6C(2) 



In— - - -L\-- 
771^ 3 3 / a 



LJln— -- --L2_5^(2) 



,(48) 



where Lb, Lc, Lbc and J(a:, y) were defined in Ea. (|44p . In the previous formula the ie prescription is implicit 
in the squared lepton mass {mf —> m,f — ie) and gives the rule to extract the proper imaginary part of the 
logarithms. For electron loops the vertex correction differs by a constant [T4] : 



Ve{q') = - 



a 
3^ 



lln34-^ln^4 + i 
12 m2 24 m? 2 



C(2)- 



265 
"36" 



1 -r 3 ,,„, 383 
ln^ + -C(2)- — 

7TT.J 4 36 



(49) 
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the remaining corrections are identical. 

In order to obtain the total leptonic corrections, the contributions of the one- and two-loop vacuum 
polarization function have to be added: 



dai daV daf 



dn dn dn dn 



d(T\°^ 



(50) 



The second term can be obtained from Eq. ([2]) with the substitution 11 -^ H; , while last term can be computed 
taking from the literature the expression for the leptonic contribution to the two-loop vacuum polarization 
function (Fig. [5]): 



daf o? ri-2a; + 2x2 



'I 

dn 



2 Renf ^ (s) 



2-2a; + a;^ 
2~^ 



2U\^\t) 



\-2x + x^ 



RenJ'^(s)+nj'^(s) 



where in the high-energy limil|j 



nf^(g^) 



a 

4^ 



In 



mf — ie 6 



4C(3) 



(51) 



Comparing our analytical result with 0|, we find perfect agreement. 



^-c 2Re 




Figure 6: Contributions involving the two-loop vacuum polarization from a lepton I. 



5 Numerical analysis 

To arrive at a numerical result we adopt the following parametrizations for R(s): For the comparison 
with earlier work [9] we take the function provided by H.Burkhardt [10]. This parametrization (denoted by 
B) is simple and efficient for the integration, however, it includes only data more than 20 years old. A newer 
parametrization (denoted by HMNT) is based on the most recent and accurate data and will be used for 
most of our detailed predictions. The two parametrizations for R{s) are shown in Fig. [T] 

The contributions from narrow resonances are incorporated using: 



Ri-csis) 



9tt 



a^iM.c 



-M,,sT,+ ,-S{s-Ml,). 



(52) 



For the parametrization HMNT we take J/^, ■0(25), T(15), T(25) and T{3S) as narrow resonances with 
the parameters listed in Table [1] thus replacing their rapidly varying cross section governed by a narrow 
Breit-Wigner shape with an easy to be integrated delta function. Parametrization B uses slightly different 



M(GeV) 
ree(keV) 

(a/a(M))2 



J/* 



3.096916(11) 
5.55(14) 
0.957785 



*(2S') 



3.686093(34) 
2.48(6) 
0.95554 



T(IS') 



9.46030(26) 
1.340(18) 
0.932069 



T{2S) 



10.02326(31) 

0.612(11) 

0.93099 



T(3S') 



10.3552(5) 
0.443(8) 
0.930811 



Table 1: Masses and electronic widths of the narrow resonances and effective electromagnetic coupling at the appro- 
priate scales. 



^For general mf/\q'^\ the result can be found in |15| 
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Figure 7: The two parametrization B (dashed) and HMNT (solid) for R{s) without narrow resonances. 



values and includes in addition w(782), $(1020), V'lSTTO), V(4040), V(4160), V(4415), T(45'), T(10860) and 
T(11020) as narrow resonances and we adopt the parameter values listed in the code [TOl. For later use we 
also give the results for the moments i?(oo), Rq, Ri and i?2 based on parametrization B: 



i?(oo) = 4.0, 



Rc] — — 8.31, 



i?i = 13.1, 



i?2 = -15.6. 



(53) 
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parametrizations B (dashed) and HMNT (solid) and the high energy approximation (dotted) 
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The results for the vacuum polarization Tl{q^) and the vertex correction V{q'^) for space-like and time-like 
q^ are shown in Fig.[S]as functions of q^. We display the predictions based on both parametrization B (dashed) 
and HMNT (solid). For comparison we also show the behaviour in the high energy approximation (dotted) 
of eg . (|inillT|) . for parametrization B only. As expected from the comparison in Fig. El the difference between 
the two parametrization leads to differences in n(g^) and V{q'^) of less than 10% which are unimportant 
for the two- loop analysis (the present uncertainty for HMNT amounts to typically one to two percent). 
For n(t) and V{q'^) the high energy approximation starts to deviate significantly from the full result for 
energies below 3 GeV, while for n(s) the resonant behaviour cannot be reproduced by this approximation. 
The corresponding results for the functions B(a,b,c), which govern the behaviour of the irreducible box 
contribution are shown in Fig. [9] for a set of representative energies as functions of the scattering angle 9. 
The result for B{s,u, t) is obtained from B{s,t, u) through the substitution cos f? — > — cos 6. 




Figure 9: The functions B(a,b,c) defined in Eq. 
tions B (dashed) and HMNT (solid). 



for different kinematical regions versus cos 9 using the parametriza- 
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As expected, the predictions for 11, V and B based on the two paratrizations B and HMNT are always 
quite close, hence the following discussion will be based on HMNT only. 

The corrections for the differential distributions are shown in Fig. [TO] for four characteristic energies, 
normalized relative to the Born predictioro. They are separated into those from reducible diagrams (dcr'^'"^ = 
da^^ + da^^), irreducible vertex (dcr^) and box (da^) diagram^. In most of the cases the reducible ones 
are significantly larger than the irreducible ones, a consequence of their enhancement by the large logarithm 
ln(s/TOg). 
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Figure 10: Relative corrections to the cross section from irreducible boxes dcr^ (solid), vertices da^ (dashed) and 
reducible contribution da'^'^'^ = dof^^ + daY^d (dotted) for four characteristic energies using parametrization HMNT. 

In Fig. [TT] we display the corresponding contributions from muons (solid line) and r leptons (dashed 
line), which can be evaluated similar to the hadronic ones. It is interesting to observe that the high energy 
approximation (hea) for the muon case, rni <^ s, |i|, |m|, (dotted) fails quite badly for small angles at ^/s — 1 
GeV, a fact that could be anticipated already from Fig. [H which shows the pour quality of this approximation 
for -s/^ < 3 GeV in the hadronic case. For high energies, the quality of the approximation should be sufficient 
for all practical purposes (Fig. fTTj) . 



•^We do not present the two-loop vacuum polarization insertions of Eq.Q wiiicli are best combined witii tiic one- loop and 
Born contribution in the resummed form of Eq. (|6ll. It is clear that nhad(9'^) niust be known with a relative precision of about 
one percent, if one aims at luminosity determination with an error significantly below one per mille. 

■^Here and below the infrared-sensitive contributions proportional to ln(2a;/y^) are set to zero. 
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Figure 11: Leptonic cross section dai from muons (upper) and r-leptons (lower) (Ea. (|45|l l relative to the Born cross 
section. Comparison is shown between the high energy approximation (hea) of Ea. (|48p (dashed) and the exact result 
(sohd) obtained from numerical integration using Ri as defined in Eq. (|46|) . 



The relative corrections from hadron and lepton (e, /i, r) loops are compared in Fig. [T^ A markedly 
different energy and angular dependence is observed for the four contributions. Individually and in the sum, 
they significantly exceed the level of one per mille necessary to achieve the corresponding precision of the 
luminosity measurements. However, as discussed before, the reducible terms dominate and the irreducible 
hadronic terms are typically below one per mille. For precise comparisons the numerical results are also 
listed in Table [H for a selected set of energies and angles. For small angles the box contribution dcP remains 
tiny, often around or below 10~^ of the Born cross section, and the result is dominated by the reducible 
correction da^^^ which is typically a factor 10 to 100 larger and is trivially obtained from existing one-loop 
results, Eq.([T3 | -(fT8 l) . 



To illustrate the relative importance of reducible and irreducible contributions, the results for the irre- 
ducible box da^ and the sum da^^''^^ = da^ + a'^"^ are listed in Table [3l The relative contribution of da^ 
is evidently tiny. The results for da^^'^^'^ are also compared with those from [^ . For the hadronic case they 
are in good agreement, although sometimes deviating in the last of the digits listed in [S]. For the leptons 
perfect agreement is observed. 
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Figure 12: Relative corrections from hadron and lepton loops as functions of cos Q for four characteristic energies. The 
high energy approximation of Ea. (|48|) is used for electron contribution at all energies and for muons for ^fs > 10 GeV. 



In general, the corrections exhibit a fairly smooth energy dependence. However, the situation changes for 
energies close to narrow resonances. This is exemplified in Fig. [13] for two cases: around the $ and around 
the J/^ resonances for three fixed angles: 3°, 90° and 177°. The interference of the continuum amplitude 
with a Breit-Wigner enhanced correction is clearly visible. At 3° (177°), irreducible box and the reducible 
corrections are of comparable size and opposite (equal) sign, at 90° the reducible ones dominate. The vertex 
corrections are always small. Formally for the case of the J/^, treated as narrow resonance, the correction 
even diverges, and it is still extremely large if the natural width of J/^ is introduced. In practice, however, 
the cross section has to folded with the energy spread of order MeV. In this case the singular amplitude with 
its asymmetric behaviour around ^/s = Mj-es is damped and thus remains a small correction. From these 
considerations, it is clear that a precise parametrization of R{s) is required in regions of rapidly varying cross 
section, if one aims at a precise prediction of the corrections in this region. 
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0.001579 
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0.03537(1) 


-3.418 
-3.200 


0.002188 
0.002055 


-0.1682 
-0.1557 


-7.017 
-6.589 


[10-3] 


-0.08749 
-0.08347 


-1.458 
-1.375 


-1.516 
-1.430 


-0.4614 
-0.4373 


-3.715 
-3.495 


-3.970 
-3.734 


da'^^^/da" 
[10-3] 


1.650 
1.562 


8.340 
7.835 


8.001 
7.526 


4.594 
4.289 


13.35 
12.60 


12.84 
12.12 


dOhad/rfcr" 
[10-3] 


1.565 
1.480 


6.918 
6.495 


3.066 
2.895 


4.135 
3.854 


9.469 
8.944 


1.855 
1.795 


dae/dcr° 
[10-3] 


0.8128 


-0.4993 


-6.769 


0.5561 


-1.954 


-10.33 


dcTp/dcr" 

10-3] 


0.7078 


1.943 


0.3714 


1.378 


2.366 


-0.3557 


[10-3] 


0.04932 


1.246 


0.9548 


0.4787 


2.100 


1.162 



Table 2: Born cross section and relative corrections from hadronic irreducible boxes and vertices, reducible hadronic 
contributions and their sum, for selected energies and angles. Upper/lower lines: parametrization B/HMNT. Also 
shown are the corresponding contributions from electrons, muons and r-leptons. 
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da^/dn 


12.8008 
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-0.00115553 


-0.000171183 


M 


da^+""^/dfl 


160.197 
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6.08187 
6.08 


0.147046 
0.1470 


0.00725789 
0.00726 




da^/dn 


0.465857 


0.0939460 


0.00188681 


0.0000195543 


r 


daB+rod/dn 


2.38272 
2 


1.33347 
1.33 


0.0752669 
0.0752 


0.00457124 
0.00457 



Table 3: Corrections from the irreducible boxes (first line), sum of box and reducible contributions (second line) and 
comparison with 9 (third line), for hadrons, muons and r-leptons at selcted energies and angles. In the upper table 
the numbers are in units of 10 nbarn, in the lower one in units of 10~ nbarn. 



The implementation of these results in a Monte Carlo generator is straightforward and their modular 
structure should lead to an efficient program. The reducible contribution da'^°'^ can be obtained from the 
one-loop corrections simply modifying the photon propagators outside the loop according to: 



I^^[l + Ren(g2)], 



n = n 



had 



He + n, 



n., 



(54) 



and adding the terms proportional to Imn(s) multiplied by the imaginary part of the one-loop result. The 
irreducible vertex corrections V{q'^) can be directly combined with the Born cross section. All these are 
one-dimensional functions that can be tabulated once for ever. The irreducible box contribution is decom- 
posedinto terms proportional to n(s) and n(i) plus a remainder characterized by the functions B(s,Lu), 
B{s, u, t), B{t, s, u) and B{t, u, s). These are obtained through efficient and precise integration routine^ 



^available upon request from the authors. 
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Figure 13: Behaviour of the hadronic corrections around the $ and the J/'^ resonances for three characteristic angles. 
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6 Conclusions 

Using published one- loop results, a compact formula has been derived which, in combination with dis- 
persion relations and the by now well-measured R-ratio, can be used to evaluate the hadronic contributions 
to Bhabha scattering. The same approach is applicable for leptonic contributions, in particular from muons 
and T-leptons. The method and result are valid in the limit m^ <C s, |t|, |u| for arbitrary R{s) and arbitrary 
m^ ^/s. Comparing with [9], our numerical results are in perfect agreement for massive leptons, with mf 
arbitrary, while for hadronic contributions we observe small numerical differences. In the high energy limit 
the integrals can be evaluated in analytic form and the results have been compared with those for lepton 
loops that can be found in the literature. We find that overall size of the corrections, their sign and their 
angular dependence differ significantly between hadron, muon, r-lepton and electron contributions. The size 
of the hadronic corrections varies from a fractional up to several permille. However, these are dominated 
by the reducible ones, with the irreducible box and vertex terms being typically below one permille. The 
modular structure of the results allows for a simple implementation into any Monte Carlo generator. For 
such an implementation, the corrections from virtual plus soft real photon radiation must be complemented 
by hard real radiation. This part is evidently straightforward, since it involves tree-level diagrams only, with 
the photon propagator dressed by hadronic vacuum polarization. 
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